Half-filled 2D bilayers in a strong magnetic field: Revisiting the v — 1/2 fractional 

quantum Hall effect 
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We examine the quantum phase diagram of the fractional quantum Hall effect in the lowest 
Landau level in half-hlled bilayer structures as a function of tunneling strength and layer separation. 
Using numerical exact diagonalization we investigate the important question of whether this system 
supports a fractional quantum Hall effect described by the non-Abelian Moore-Read Pfafhan state 
in the strong tunneling regime. We find that, although it is in principle possible, it is unlikely that 
the non-Abelian FQHE exists in the lowest Landau level. We establish that all so far observed 
FQHE states in half-filled lowest Landau level bilayers are most likely described by the Abelian 
Halperin 331 state. 
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Two important developments have rekindled interest 
in the phenomena of even-denominator incompressible 
fractional quantum Hall (FQH) states in bilayer semi- 
conductor structures. The first is the recent intriguing 
experimental observation by Luhman et al. [l[ of two 
distinct even-denominator (filling factors v = 1/2 and 
1/4) FQH states in a very wide (~ 600A) single quan- 
tum well (WQW) structure at very high (> 40T) mag- 
netic fields. The second development, motivated by im- 
plications for fault-tolerant topological quantum compu- 
tation [2J, is the great deal of recent theoretical and ex- 
perimental interest in the possible non-Abelian nature of 
the v = 5/2 second Landau level (SLL) FQH state in 
the single-layer system This leads to the interest- 
ing as well as important question whether a non-Abelian 
v = 1/2 FQH state, i.e., the analog of the possibly non- 
Abelian v = 5/2 = 2 + 1/2 SLL single-layer state, can 
exist in the lowest Landau level (LLL) under experimen- 
tally observable conditions. This question has a long his- 
tory 0] in the theoretical literature going back to the 
early 1990s. Motivated by these considerations we revisit 
the v = 1/2 FQHE in bilayer structures (with each layer 
having v = 1/4 on average) by theoretically investigating 
the quantum phase diagram for the LLL v = 1/2 bilayer 
FQH state using the FQH spherical system finite size ex- 
act diagonalization technique. We obtain an approximate 
quantum phase diagram for the spin-polarized v = 1/2 
bilayer FQH system in the LLL in the parameter space of 
inter-layer separation (d), inter-layer tunneling strength 
(t) or, equivalently, the symmetric-antisymmetric gap en- 
ergy, and the width or thickness (w) of the individual 
wells (for simplicity, we assume the two wells to be iden- 
tical). We take the system to be fully spin-polarized since 
all experimental v = 1/2 FQH states appear to be fully 
spin-polarized. 



Our main findings are 
WQW v = 1/2 FQHEs 0. 
Halperin 331 FQH states 



(i) the recently observed 
are strong-pairing Abelian 
which, however, sit close to 



the boundary between the Abelian 331 and the weak- 
pairing non- Abelian Pfaffian (Pf) FQH state 0, and (ii) 
it may be conceivable, as a matter of principle, to real- 
ize the LLL v =1/2 Pfaffian non-Abelian FQHE in very 
thick bilayers, but as a matter of practice, this is unlikely 
since the v =1/2 FQHE gap is extremely small (perhaps 
zero) in the parameter regime where the Pf is more sta- 
ble than the 331 phase. Our findings about the fragility 
of the LLL v = 1/2 non-Abelian Pf state are consistent 
with recent conclusions 0, Q , but our main focus in the 
current work is in understanding the v = 1/2 bilayer 
quantum phase diagram treating t, d, and w as indepen- 
dent tuning parameters of the system Hamiltonian. 

Before presenting our results, we discuss the context 
of our theoretical investigation. The proposed Moore- 
Read Pf state is a weak-pairing single-layer FQH state 
at v = 1/2 which, in principle, applies to any orbital LL 
(i.e., LLL as well as SLL). Thus, as a matter of principle 
the v =1/2 LLL single-layer Pf FQHE is certainly a pos- 
sibility although it has never been observed experimen- 
tally. The best existing numerical work 0, [§] indicates 
that either the single-layer v = 1/2 LLL Pf state does 
not exist in nature or if it exists, does so only in rather 
thick 2D layers with an extremely small FQH excitation 
gap, making it impossible or very difficult to observe ex- 
perimentally. By contrast, the single-layer v = 5/2 SLL 
FQHE is observed routinely, albeit at low temperatures 
(< lOOmK), in high mobility (> 10 7 cm 2 /Vs) samples, 
and with a rather small (but experimentally accessible) 
activation gap (~ 100-500mi4T). In fact, it has been 
pointed out that the experimental v = 5/2 FQHE is al- 
ways among the strongest (along with the v = 7/3 and 
8/3 FQHE) observed FQH states in the SLL. 

Instead of studying a single-layer 2D system we con- 
centrate on the spin-polarized bilayer system assuming 
an arbitrary tunneling strength t (proportional to the 
symmetric-antisymmetric gap) and an arbitrary layer 
separation d. We numerically obtain the quantum phase 
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diagram at v = 1/2 for this bilayer system in the t-d 
space, concentrating entirely on the Pfaffian and the 331 
FQHE phases. Our focus on the bilayer spin-polarized 
v = 1/2 FQHE is consistent with the experimental fact 
that the incompressible v = 1/2 FQHE has so far been 
observed only in effective bilayer structures. The goal 
is to carry out an extensive comparison with all existing 
bilayer v = 1/2 FQH experimental observations to ascer- 
tain any hint of a Pfaffian v = 1/2 state for large values 
of t. It is obvious that our model system is an effective 
bilayer (single-layer) v = 1/2 system for small (large) val- 
ues of t, and therefore by studying the quantum phase 
diagram as a function of t (and d) we hope to shed light 
on the possible existence of a single-layer v = 1/2 FQHE 
in real systems. Recent FQHE experiments at v = 1/2 
for large values of the tunneling strength make it imper- 
ative that a theoretical analysis be carried out to achieve 
a proper qualitative understanding of the experimental 
situation [3, 0|. 

We use a simple model Hamiltonian incorporating both 
finite tunneling and finite layer separation for our bilayer 
FQH system: 

N 

H = ^[MntraCI^-rj-D+^tradri-fjl) 

i<j 

+ VWIr, . (1) 

where r^(?i) is the position of the i-th electron in the 
right (left) layer. In E q. Q, V i ntaa (r) = e 2 /(«Vr 2 + w 2 ) 
and VintcrM = e 2 /{n\fr 2 + d 2 ) are the intralayer and in- 
terlayer Coulomb interaction incorporating a finite layer 
width w and a center-to-center interlayer separation d 
(> w by definition). The x-component of the pseudo- 
spin operator S x controls the tunneling between the two 
quantum wells with large t denoting strong tunneling. 
We numerically diagonalize H in the spherical geome- 
try assuming specific values of w, d, and t (each ex- 
pressed throughout in dimensionless units using the mag- 
netic length I = {cfi/ eB) 1 / 2 as the length unit and the 
Coulomb energy e 2 /(nl), where k is the background di- 
electric constant, as the energy unit). Following the stan- 
dard well-tested procedures [J, d, Q used extensively in 
the FQHE literature, we calculate the overlap between 
the exact numerical N electron ground state wavefunc- 
tion of the Coulomb Hamiltonian defined by Eq. |T]) and 
the candidate N electron variational states which are the 
Halperin 331 strong-pairing Q and the Moore-Read Pfaf- 
fian weak-pairing [7| wavefunctions: 

N/2 N/2 N 

*331 = liiZi-Z^lliZi-ZjflliZi-Zj) (2) 

i<j i<j i,j 

* Pf = Pf{^:}nfe-^) 2 . (3) 

^ 1 3 ' i<j 

respectively, where z = x— iy is the electron coordinate in 



the plane. We ensure that (i) the ground state is homo- 
geneous, i.e., has zero total orbital angular momentum, 
and (ii) there is a gap, the FQH excitation gap, separat- 
ing the ground state from all excited states. The results 
shown in this paper all use an N = 8 electron system, 
but larger systems show the same qualitative features. 
Since the theoretical techniques are standard, we do not 
give the details, concentrating instead on the results and 
their implications for v = 1/2 FQHE experiments. 

The calculated overlap and gap determine the nature 
of the FQHE and its strength at v = 1/2 in our the- 
ory. We operationally define the system to be in the 331 
(Pf) phase depending on the overlap with the 331 (Pf) 
state being the larger of the two. We emphasize that our 
work is a comparison between these two incompressible 
states only, and we cannot comment on the possibility 
of some other state (i.e., neither 331 nor Pf) being the 
ground state. We do, however, believe that if the system 
is incompressible at a particular set of parameter values 
(i.e., d, t, w, etc.), it is very likely to be described by one 
of these two candidate states, 331 or Pf, since no other 
proposed incompressible state exists in the literature for 
v = 1/2 spin-polarized bilayer or single-layer LLL sys- 
tems. We cannot, however, rule out the possibility that 
the real system has a compressible ground state (with- 
out manifesting FQHE), e.g., a composite fermion sea, 
not considered in our calculation. 

We first show in Fig. [1] our numerically calculated 
FQHE quantum phase diagram (QPD) for the bilayer 
v = 1/2 system in the t-d space with the color coding in- 
dicating the numerical FQH gap strength and the dashed 
line separating the 331 phase from the Pf phase (i.e., 
the overlap with 331 (Pf) being larger above (below) the 
dashed line). Note that the dashed line is only an oper- 
ational phase boundary within our calculation since all 
we know is that the 331 (Pf) has higher (lower) overlap 
above (below) this line. The calculated overlap for each 
phase, i.e., 331 (Pf) above (below) the dashed line, varies 
between ~ 0.6 and ~ 0.96 for our 8-particle system. We 
show the QPD for three values of the layer width pa- 
rameter (a) w = 0; (b) w = 0.6; (c) w = 2.4. The 
zero- (Fig. [TJa)) and the intermediate-width (Fig. |Tfb)) 
results are of physical relevance whereas the (unrealisti- 
cally) large width results (Fig. [T^c)) are provided here 
only for completeness (since this is the regime where the 
Pf state dominates over the 331 state in the QDP). We 
note that we are using the simplistic Zhang-Das Sarma 
(ZDS) model for describing the well width effect, and 
crudely speaking w = 1 in the ZDS model corresponds 
roughly to wqw ~ 6 where wqw is the corresponding 
physical quantum well width. For a single WQW, where 
the effective bilayer is created by the self-consistent po- 
tential of the electrons themselves, our model w is typi- 
cally much less than the total width W of the WQW-very 
roughly speaking w ~ W/6, and d ~ W/2. As empha- 
sized above, we treat t, d, and u>(< d) as independent 
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FIG. 1: (color online) Quantum phase diagram (QPD) and FQHE gap (color coded) versus layer separation d and tunneling 
strength t for widths (a) w — 0, (b) w = 0.6, and (c) w = 2.4. For the QPD, the 331 and Pf phases (as discussed in the text) are 
separated by a dashed black line and labeled appropriately. The FQHE gap is given as a contour plot with color coding given 
by the color-bar from dark to light, i.e., white being a largest value of 0.4 and black being value of 0. The asterisks, triangles, 
circles, and squares correspond to the different experiments in Refs. Il2l . [Til . @ and [JJ, respectively. Only experimental points 
showing FQHE are within the large solid circles in (a) with the lower smaller (upper larger) circles indicating experiments in 
double-quantum-well (single- WQW) structures. We note that the single triangle on the Pf side of the QPD does not manifest 
any experimental FQHE indicating that the theoretical gap may be overestimated for the Pf state. It is obvious from the 
figures that the Pf state becomes more dominant, albeit with very small FQH excitation gap, for larger values of w. 



tuning parameters. 

In Fig. [TJ we have put as discrete symbols all exist- 
ing v — 1/2 bilayer experimental data (both for double 
quantum well systems and single wide quantum wells) in 
the literature, extracting the relevant parameter values 
(i.e., d and t) from the experimental works [U IH, [ill \Y2\ . 
Because of the ambiguity and uncertainty in the defini- 
tion of w, we have put the data points on all three QPDs 
shown in Fig. [TJ although the actual experimental width 
values correspond to only Figs.QJa) an d (b). 

Results shown in Fig. [TJ bring out several important 
points of physics not clearly appreciated earlier in spite of 
a great deal of theoretical exact diagonalization work on 
v = 1/2 bilayer FQHE: (i) It is obvious that large (small) 
t and small (large) d, in general, lead to a decisive pref- 
erence for the existence of v = 1/2 Pf (331) FQHE. The 
fact that large t values would preferentially lead to the 
Pf state over the 331 state is, of course, expected since 
the system becomes an effective one-component system 
for large tunneling strength, (ii) What is, however, not 
obvious, but apparent from the QPDs shown in Fig.[TJ is 
that the FQH gap (given in color coding in the figures) 
is maximum near the phase boundary between 331 and 
Pf. (iii) Another non-obvious result is the persistence of 
the 331 state for very large (essentially arbitrarily large!) 
values of the tunneling strength t as long as the layer 
separation d is also large-thus having a large t by itself, 
as achieved in the Luhman et al. experiments [l[, is not 
enough to realize the single-layer v — 1/2 Pf FQHE, one 
must also have a relatively small value of layer separa- 
tion d so that one is below the phase boundary (dashed 
line) in Fig. [TJ The explanation for the Luhman experi- 
mental v = 1/2 FQHE being a 331 sate, as can be seen 
in Fig. [TJ is indeed the fact that both t and d are large 



in these samples making 331 a good variational state, 
(iv) An important aspect of Fig.[TJis that the Pf FQHE 
gap tends to be very small-this is particularly true for 
larger values of w, where the Pf overlap is larger. This 
implies, as emphasized by Storni et al. [9J, that the obser- 
vation of a v = 1/2 Pf state is unlikely since the activa- 
tion gap would be extremely (perhaps even vanishingly) 
small, (v) For larger values of w (and large t), our cal- 
culated QPD is dominated by the Pf state-particularly 
for the unrealistically large width w — 2.4 (correspond- 
ing to wqw ~ 14!) where all the experimental d and 
t values fall in the Pf regime of the phase diagram. We 
emphasize, however, that this Pf-dominated large-w (and 
large-i) regime will be difficult (perhaps even impossible) 
to access experimentally since the FQH gap would be ap- 
parently extremely small as in Fig. [TJ 

In discussing Fig. [TJ further, we mention that our 331 
(Pf ) regimes not only have the wavefunction overlap with 
the corresponding 331 (Pf) state being larger than the 
other, but also the calculated expectation value (S x ) rs 
(4) in the 331 (Pf) regime. Thus, our QPD is consistent 
with both the overlap and the pseudo-spin calculation as 
obtained from exact diagonalization. 

We now discuss the published experimental results in 
light of our theoretical QPD. First, we note that most 
of the existing experimental points fall on the 331 side 
of the phase diagram which is consistent with our QPD 
in Fig. [TJ In particular, only samples on the 331 side of 
the QDP with reasonably large FQH gaps, i.e., the data 
points close to the phase boundary, exhibit experimental 
FQHE. By contrast, the one data point (in Figs.QJa) and 
(b)) on the Pf side of the phase boundary does not man- 
ifest any observable FQHE in spite of its location being 
in a regime of reasonable FQH excitation gap according 
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FIG. 2: (color online) FQHE energy gap versus (a) tunneling 
strength t for a few values of layer separation d and (b) layer 
separation d for a few values of t (both consider zero width 
w = 0) . A dashed vertical line of the same shade corresponds 
to the boundary between the PfafHan phase (right/left of the 
line for (a)/(b)) and the 331 phase (left/right of the line for 
(a)/(b)). 



to our phase diagram. This is consistent with the find- 
ing of Storni et al. [9] that the v = 1/2 FQH Pf gap in 
a single-layer system is likely to be vanishingly small in 
the thermodynamic limit. It is, therefore, possible that 
much of the Pf regime in our QPD has a much smaller 
excitation gap than what we obtain on the basis of our 
N = 8 particle diagonalization calculation. We refer to 
Storni et al. § for more details on the theoretical status 
of the single-layer LLL v = 1/2 FQHE. 

For a more detailed view of the v = 1/2 bilayer FQHE, 
we show in Figs. Ufa) and (b), respectively, our calculated 
FQHE gap as a function of t (for a few fixed d values) 
and as a function of d (for a few fixed t values). In each 
figure, we also depict the line separating the 331 (smaller 
i/larger d) and the Pf (larger i/smaller d) regimes in the 
phase diagram. The qualitatively interesting point is, of 
course, the non-monotonicity in the FQHE gap as a func- 
tion of t or d with a maximum close (but always on the 
331 side) to the phase boundary. The non-monotonicity 
in the FQHE gap as a function of t (but not d) was earlier 
pointed out, but our finding that the peak lies always on 
the 331 side of the phase boundary is a new result. We 
emphasize that the FQHE gap peak lying always on the 
331 side of the phase boundary is strong evidence that 
the 331 phase is the dominant FQH phase in v = 1/2 



systems. We believe that the only chance of observing 
the v = 1/2 Pf FQHE is to look on the Pf side of phase 
boundary at fairly large values of d and t. This is in sharp 
contrast to the SLL v = 5/2 bilayer FQHE where we re- 
cently showed that there are two sharp ridges far away 
from each other in the d-t space corresponding to the 
v = 5/2 Pf and 331 bilayer phases fl3| . We note that for 
unrealistically large w (Fig. [He)), Pf dominates over 331 
but the FQHE gap becomes extremely small everywhere. 

We conclude by commenting on the nature of the quan- 
tum phase transition (QPT) between the 331 and the 
Pf phase in the d-t space. It may appear at first sight 
that our work implies a continuous QPT from the strong- 
pairing 331 to the weak-pairing Pf state with increasing 
(decreasing) t (<f). This is, however, deceptive since all 
we are doing is comparing these two phases using a fi- 
nite size diagonalization study for discrete values of t 
and d. It is entirely possible that a completely different 
phase, e.g. a compressible composite fermion Fermi liq- 
uid phase, has lower energy and intervenes between the 
331 and Pf phases so that the system goes from 331 to Pf 
(or vice versa) through two first-order transitions. There 
is independent numerical evidence that the compressible 
composite fermion sea indeed has a lower ground state 
energy than the Pf state in a single-layer v = 1/2 (but 
not 5/2) system which corresponds to the large t (and 
small d) regime in our QPD. This would indicate first 
order transitions in going from 331 to the Pf (if it exists) 
through the compressible phase. What we have shown 
here is that if the v = 1/2 bilayer Pf phase exists at all, 
then it would manifest most strongly in wide samples and 
close to the phase boundary with the 331 phase, but will 
have an extremely small FQHE excitation gap. 

This work is supported by Microsoft Q and DARPA 
QuEST. 
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